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Abstract: The paper highlights the investigation of nonlinear axial-force 
and moment-curvature relationships for rectangular reinforced concrete 
(RC) beam and column cross sections for the reference reinforced concrete 
building structure considered for Probabilistic Seismic Risk Evaluation. The 
highly popular model namely Mander’s model is used for concrete stress- 
strain relationship since it is simple and effective in considering the effects of 
confinement. The module determines the expected behavior of a user-defined 
cross-section by first dividing the section into a number of parallel concrete 
and steel “fibers”. Then, the section forces and deformations are determined 
from the fiber sfrains and stresses using fundamental principles of equilibrium, 
strain compatibility, and constitutive relationships assuming that plane sections 
remain plane. 

Index Terms: Mander’s Method, Moment-Curvature, Stress-Strain Relationship 
and Moment Curvature 


1. INTRODUCTION 


Since the country lie in earthquake prone area and many of the destructive 

earthquakes occurred in the history so far resulting in high number of 

casualties due to collapse of buildings and dwellings. A major challenge for 

the performance based seismic engineering is to develop simple yet efficiently 

accurate methods for analyzing designed structures and evaluating existing Journal on Today’s ideas- 

buildings to meet the selected performance objectives Elastic analyses are Tomorrow s'^chnoiogies, 

insufficient because they cannot realistically predict the force and deformation 2014 

distributions after the initiation of damage in the building. Inelastic analytical pp. 13-29 
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procedures become necessary to identify the modes of failure and the potential 
for progressive collapse. The need to perform some form of inelastic analysis 
is already incorporated in many building codes. Theoretical moment-curvature 
analysis for reinforced concrete columns, indicating the available flexural 
strength and ductility, can be conducted providing the stress-strain relation 
for the concrete and steel are known. The moments and curvatures associated 
with increasing flexural deformations of the column may be computed for 
various column axial loads by incrementing the curvature and satisfying the 
requirements of strain compatibility and equilibrium of forces. 


2. STRESS-STRAIN MODEL EOR CONCRETE 

Mander’s model is highly popular model since it is simple and effective in 
considering the effects of confinement. It considers increase in both the strength 
and ductility of RC members with confined concrete. The model is popularly 
used to evaluate the effective strength of the columns confined by stirrups, 
steel jacket and even by FRP wrapping as accomplished in Figure 1. 



Figure 1: Mander’s Model for Stress-Strain Relationship for Confined 
Concrete 


3. STRESS-STRAIN CURVES FOR REINFORCING STEEL 


The idealized stress-strain curve for concrete as recommended by IS: 456-2000 
is as shown in Figure 2. Stress-strain curve for steel as per British code CP 110- 
1972 as shown in Figure 3, accordingly the term 0.7fy is the simplification of 


the expression 


fy 


7 + 

' m ' 


fy 

2000 . 


. It gives ah the simplified general equations which 


can be used for any grade of steel. 
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Figure 2: Stress Strain Curve for Concrete 


Moment Curvature 
Characteristics for 
Structural Elements 
of RC Building 



Figure 3: Stress-Strain Curve for Steel 


4. STRUCTURAL SYSTEM 

The building is an RC G-t-3 framed structure. The floor plan is same for all 
floors. The beam arrangement is different for the roof. It is symmetric in both 
the direction. The concrete slab is 120 mm thick at each floor level. Overall 
geometry of the structure including the beam layout of all the floors is as shown 
in Figure 4. 

Figure 5 and 6 shows the size and reinforcement details for floor beam 
and roof beam sections at the column face. Figure 7 shows the size and 
reinforcement details for column at the beam face. 

5. MOMENT CURVATURE EOR BEAM 

The most fundamental requirement in predicting the Moment Curvature 
behaviour of a flexural member is the knowledge of the behaviour of its 
constituents. With the increasing use of higher-grade concretes, the ductility 
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Figure 4: Overall Geometry of the Structure 
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Figure 5: Details of Roof Beams 


of which is significantly less than normal concrete, it is essential to confine 
the concrete. In a flexure member the shear reinforcement also confines the 
concrete in the compression zone. The relationship for the bending member as 
depicted in Figure 8 is as follows, 

M _ _ 1 _ 

El E a: R 
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Figure 6: Details of Floor Beams 
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Figure 7: Details of columns at various levels 
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Where, 

= Curvature 

fc = Maximum stress in compression in concrete 
e = Maximum strain in concrete 

C 

e = Maximum strain in steel 

S 

X = Depth of neutral axis 
d = Effective depth of section 

Let consider BF205 beam section as depicted in Figure 9, 




Figure 9: Detailed BF204 Beam Section 


280 

Modular Ratio (m) for elastic analysis = -= 14 

4 

Curvature at cracking moment just before cracking 


f = 0.7,/f7 = 3.1305 MPa 

At just before cracking moment is resisted by concrete in tension therefore 
Neglecting steel, N.A. = d/2 

f 1 

M = = 120.00 KNm 

cr 

y 

f M 

ip = ^ = — = 0.00028 

y El 

Curvature after cracking 
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Concrete does not take any tension so depth of neutral axis x is, 


k = ^\(,p +p'fn +2n{p +p'‘^)-{p +p')n 


k = 0.2179 


Effective moment of inertia 


bx^ 

I „ =-+ A m(D - x)" = 7.990mm'‘ 
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$atM =—2^ = 0.00067 

“ EIe« 

£ Strain in steel = f /E = 0.0021 

s y s 

The neutral axis at yielding is given as distance kd from extreme 
compression fiber, where the ratio k is calculated using expression: 


Mp + p'f n +2n{p + p'^)-{p + p')n 


k = 0.2179 


Where p 


and p' 


A, 

s 

bd 


are the tension and compression steel ratios, n is 

the modular ratio, and d and d’ are the distance of compression and tension 
steel from extreme compression fiber. 

A Steel in tension 
A^' Steel in compression 

Depth of Neutral axis at yield stage; x = kd =208.65mm = .21m 

Taking moment about compressive force due to concrete, yield moment 
is given by: 


M =AJ 


kd 

d - 

3 




I r! kd , 

/- d 

s •’ y 


Since stress in the tension steel is f^, using similar triangles, stress in compression 
steel is calculated as 

7 if 

f = = 99.01MPa 

d-kd ^ 

Curvature is then obtained as 

If = = 0.00277 

“ d - kd 
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After yielding of tension steel, its stress remains constant but strain keeps 
increasing until compressive strain in extreme fiber of concrete reaches tbe 
strain value of e at maximum stress in concrete f ’. In order to address tbe 

cu c 

nonlinearity in concrete at bigb strains, Wbitney-block is used to approximate 
tbe parabolic stress distribution in concrete to an equivalent rectangular stress- 
block representation. 

The calculation of ultimate state requires iteration. For hand calculations, 
let us assume that strain in compression steel e^, exceeds the yield strain e^. 
This assumption will be checked later. Equilibrium of tension and compressive 
force gives the depth of neutral axis as; 


A,f,-AX 

85f>/3, 


A f — a' f' 

——^ = 187.3319 
36f\b 


Ultimate moment is then obtained by taking moment about tension steel as: 

c 


M = C 


d-- 

2 


= 36/' Cb 
= 361.57 KNm 


+ C^id-d') 

+ A'j\{d-d') 


d-- 

2 


The ductility factor 


^ 00^^0 01868 


D. F. = ^ = 6.74 




Where is cj)^ ultimate strain in concrete at maximum stress, which is 0.0034 
as per IS456-2000. This is first trial value of e^. Assumption of yielding in 
compression is now checked by ensuring: 

T-d 

Es =(-)ecu >ey 

c 

If the above condition is satisfied then assumption made is true and obtained 
value of(M^, cj)^) defines the ultimate state on the moment-curvature curve. If 
the condition is not satisfied further iteration is required with new trial strain 
value as 


6 = 


e -I- e 

y s 


6 . CONCRETE PROPERTIES 

Currently, the entire cross-section is assumed to be unconfined. The compression 
stress-strain relationship of the unconfined concrete is determined using a 
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method developed by Mander et al. In this method, the concrete stress, fc is 
given as a function of the strain, as: 

= Jc^in MPa 
r-l + x^ 



e^gis the strain at peak stress ( f^’) and is the tangent modulus of elasticity 
of concrete calculated as 

Ec = 5000 ^/F MPa 

Esec, the secant modulus of elasticity, is the slope of the line connecting the 

f' 

origin and peak stress on the compressive stress-strain curve (i.e., Esec=-^). 

f 

CO 

Crushing of the unconfined concrete is assumed to occur at e 2e . 

^ CU CO 


7. THE MOMENT-CURVATURE RELATIONSHIPS OE THE 
COLUMN CROSS-SECTIONS 

This is an iterative process, in which the basic equilibrium requirement P = C-T 
is used to find the neutral axis location, c, for a particular maximum concrete 
compressive strain, where P = axial force; C = internal compression stress 
resultant; and T = internal tension stress resultant. The total concrete compressive 
stress resultant and the location of its centroid are determined by integrating 
numerically under the concrete stress distribution. The bending moment is 
assumed to act such that the top surface of the cross-section is in compression. 

The entire process can be summarized for a cross-section with two layers 
of reinforcing bars as depicted in Eigure 10 follows 



Eigure 10: The Cross-Section with Two Eayers of Reinforcing Bars 


Eor example the column sections property in present work is 400mm x 
900mm in size with 12-28mm dia bars placed as shown in the figure 11 below. 
The transverse reinforcement for the columns is provided 10mm stirrups/ties 
@ 100mm c/c as figure shows. 
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I. wg .1 


900 


Long. Rcinfr. 12-284> 
TnuuReinA: 10^Sl00c.‘c 



Figure 11: Detailed CL 15 and 19 Column Section (Ground to 2nd Floor) 


The calculation of the following four points on a moment-curvature curve 
will he taken for plotting the moment curvature curve for various column 
section shown in this example: e =0.25 e ;e = 0.5 e ;e = 0.75 e ; e = 
1 . 0 e^^(concrete crushing). 


A. Axial Load, P 


The axial load considered for sample calculations of CL 15 is 50% of the 
balanced failure load. The balanced load, is computed as follows: 

£ 

The neutral axis, c = = d -—, 

£ -t- £ 

cu y 

Where, e -f /E 

y y s 

With ecu = 0.004 and d = 848m, = 0.525m. 

The concrete compressive resultant, Cc, is determined by numerically 
integrating under the concrete stress distribution curve. 


cc=l fM^ = l lb 

Cc = 3531.24KN 


— d£ 


The spacing of each layer of steel is s = 0.198m 

The steel forces of each layer top to bottom wise, and F^^’ Fs 3 > ^^4 
F^j, respectively, are calculated using similarity to find the strains in the layers. 
Balanced failure condition, by definition, has strain values e = e = 0.0034 
for concrete and 6 ^^= 0.002 for bottom layer of steel. For the top steel, 

c-d' 

e=e, -= 0.0031 


Similarly for second and third layer. 
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c — d' — s 

Moment Curvature 

-= 0.0019 

// r 

Characteristics for 


Structural Elements 

^ ^ ^=0.0008 
c 

of RC Building 

d — c — 2s 


- 0.0006 



c 


—= 0.0021 
c 

This implies that the first and fifth layer of steel is at yield stress. 

Hence, 

766.61 KN (compression) 

= £s 2 ^s'^s 2 “ 458.06KN (compression) 

F^j = e^jE^A^j^lOS.TVKN (tension) 

F^^ = 142.45KN (tension) and 

E^j = fyA^ss^ 766.61KN (tension) 

Where, A^j, A^^’ A^^ are the total reinforcing steel areas in each layer. 

As per IS456 consideration the concrete tensile strength in the tension region 
recommends the modulus of rupture to he taken as, 

fr ~ 'Jj^c MPa for normal weight concrete. 

Thus, for — 20 MPa, f\ = 3-3 MPa. The concrete tension force 

C, = = 22.97 KN (tension) 

Where, A is the area of concrete in tension calculated based on the linear 
strain diagram. 

Then, the balanced axial load is found from equilibrium as, 

P =r +F +F —F —F —F —r 

'^c ^ si ^ s2 ^ s3 ^ s4 s5 '^cr 

Therefore, 50% of the balanced load used is P =1814.20kN. 

B. Instant Centroid 

The axial load acts at an “instant” centroid assumes location that for the 
calculation of the moment-curvature relationship. The location of the instant 
centroid is determined by assuming an initial condition where only the user- 
selected axial load acts on the cross-section without moment. This loading 
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condition produces a uniform compression strain distribution throughout the 
cross-section. 

Let the uniform compression strain be equal to then, 

P = fci^c + fsl^sl + ^82^52 + fs3^s3 + fs4^s4 + 

From equilibrium, 

A trial-and-error solution is needed since it is not known in advance. 
Then, the location of the instant centroid, x, from the top compression face is 
determined as, 

^ + fA.d + LAA -s) + lAjd - A + 
f A + f,A + f„A + f A + f A + f A 

•'ci c "'si .si s2 s2 ■^.s3 .s3 ''.s4 .s4 ‘'.sS s5 

X = 0.4505m 

The calculation of the first point on the moment-curvature relationship of the 
section can be summarized as follows: 


1. e =0.25e =.00084 

cm cu 

2. Assume the neutral axis depth, a distance c = .4m 

3. From the linear strain diagram geometry 

= 0.0031 

d — c 


Similarly for second and third layer, 

c — d' — s 


d-c 


= 0.0019 


e 3 = e ^ ^ = 0.0008 

c 

e =e ^ ~ ^ ~ = 0.0006 


.9 4 cm 


d — c 


= 0.0021 


4. The steel stress resultants are, 
F j = f A = 282.63 KN (comp) 


= = 84.79 KN (comp) 

F ^3 = = 129.80 KN (tension) 

Es 4 = 26.17 KN (tension) 

F^j = 350.15 KN (tension) 


24 












Determine Cc hy integrating numerically under the concrete stress 

Moment Curvature 

distrihution curve. 

Characteristics for 


Structural Elements 

Cc = £ fpdx = f bj-de^ 

of RC Building 

Cc = 1755.09 KN 



The concrete that has not cracked helow the neutral axis contrihutes to the 
tension force C^. 

/;=V7^ MPa 

= 30.82 KN 

6 . Check to see if 

Pb = + F^i + F^2 + F^3 - F^4 - F^J - C, 

P, = 1585.57 KN 

So, the neutral axis must he adjusted downward, for the particular maximum 
concrete strain that was selected in Step 1, until equilihrium is satisfied. This 
iterative process determines the correct value of c. Trying neutral axis depth c 
= .431m gives, 

= 00084 (below yield); es4 = 0.0006; = 0.0008; = 0.0019; = 0.0031 

And 

Cc = 1917.44KN; Ct = 28.51KN; 

F , =285.29KN; F , = 95.34KN; 

sl ’ s2 ’ 

F 3 =101.08KN; F ^ = 6.23KN; 

s3 ’ s4 ’ 

F =293.91KN 

s5 

Section curvature can then he found from, 

^ ^ Cm ^ 0.002125 
c 

The internal lever arms for the resultant compression and tension forces of the 
concrete measured from the instant centroid. 

Then, the section moment can he calculated as, 

M = C|,Zj. + F^jZ^j + Fj,2Z^2 + F^jZ^j + Fj^z^4 + F^^z^^ + C,Zj 

M = 872.22 kNm 

Similarly curvature and moment value are also calculated at e = 0.5e ; e = 
0.75e ;e = l.Oe (concrete crushing) 

cu’ cm cu ^ 
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The moment-curvature plot for column has been shown in Figure 12 and 
for beam as shown in Figure 13. 



Figiure 12: Moment-Curvature Relationship Curve for Column 



Figure 13: Moment-Curvature Relationship Curve for Beam 


Table 1 shows point at moment-curvature curve for column. Table 2 shows 
moment at various points in beams. Table 3 shows curvature at various points 
in beams. Table 4 shows moment at various points in column and Table 5 
shows curvature at various points in column. 
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Table 1. Point At Moment-Curvature Curve For Column 


Point 

M(moment) 

Curvature (|) 

0 

0 

0 

0.25e 

cu 

872.2165 

0.0019 

0.5e 

cu 

1395.398 

0.0046 

.75e 

cu 

1579.791 

0.0076 

l.Oe 

cu 

(concrete crushing). 

1520.344 

0.0095 


Table 2: Moment At Various Points In Beams 


Beam 

Before 

Cracking 

After 

Cracking 

At 

yield 

Crushing of 
Concrete 

BF205 

120.00 

120.00 

365.62 

361.57 

BF204 

120.00 

120.00 

153.21 

154.78 

BF225 

120.00 

120.00 

237.57 

238.17 

BF223 

120.00 

120.00 

366.72 

363.35 

BR6 

120.00 

120.00 

237.45 

238.21 

BR7 

43.20 

43.20 

88.19 

88.20 

BR21 

120.00 

120.00 

236.39 

236.97 

BR20 

120.00 

120.00 

237.57 

238.17 


Table 3: Curvature At Various Points In Beams 

Beam 

Before 

Cracking 

After 

Cracking 

At yield 

Crushing of 
Concrete 

BF205 

0.00028 

0.00067 

0.00277 

0.019 

BF204 

0.00028 

0.00138 

0.00253 

0.044 

BF225 

0.00028 

0.00096 

0.00262 

0.031 

BF223 

0.00028 

0.00067 

0.00275 

0.020 

BR6 

0.00028 

0.00096 

0.00262 

0.030 

BR7 

0.00047 

0.00159 

0.00453 

0.046 

BR21 

0.00028 

0.00096 

0.00265 

0.027 

BR20 

0.00028 

0.00096 

0.00262 

0.031 
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Table 4: Moment At Various Points In Column 


Column 

Origin 

Yield 

Ultimate 

Strain 

hardening 

CL 15/19 G/2nd 

0 

987.64 

1506.64 

802.35 

CL 15/19 3rd 

0 

440.52 

673.06 

228.49 

CL 15/19 4th 

0 

272.33 

406.20 

217.02 

CL 16/20 G/2nd 

0 

792.80 

1226.82 

647.90 

CL 16/20 3/4th 

0 

469.46 

723.75 

302.12 


Table 5: Curvature At Various Points In Column 


Column 

Origin 

Yield 

Ultimate 

Strain 

hardening 

CL 15/19 G/2nd 

0 

0.0054 

0.0078 

0.1468 

CL 15/19 3rd 

0 

0.0072 

0.0102 

0.1938 

CL 15/19 4th 

0 

0.0072 

0.0104 

0.1938 

CL 16/20 G/2nd 

0 

0.0054 

0.0073 

0.1462 

CL 16/20 3/4th 

0 

0.0054 

0.0086 

0.1281 


8. CONCLUSION 

An analytical model is presented to simulate the moment curvature behavior 
of reinforced concrete. Based on control of load increments, the algorithm 
proposed by Mander enables determination of moment - curvature-strain 
relationship with any geometry and material properties up to the maximum 
capacity of the section; however with a constant axial load or control of 
deformation increments, this model can be used to compute both the ascending 
and descending branches of the moment-curvature curve. The moment- 
curvature (or moment-rotation) relations play an important part in the study of 
limit analysis of two or three dimensional reinforced concrete frames. 


REFERENCES 

[1] Bai, Zhizhou, “Nonlinear Analysis of Reinforced Concrete Beams and Columns with Special 
Reference to Full-range and Cyclic”, Ph.D Thesis, The University of Hong Kong, China, 2006. 
http://hdl.handle.net/10722/50666 

[2] Cinitha A*, Umesha R K, Nagesh R. Iyer, “Evaluation of Seismic Performance of Existing Steel 
Buildings", American Journal of Civil and Structural Engineering, 2014, l(2):23-33. http:// 
dx.doi.org/ 10.12966/ajcse.04.02.2014 


28 











[3] Filip C. Filippou, Ahmad Issa, "Non-Linear Analysis of Reinforced Concrete Frames under 
Cyclic Load Reversals", Report No. UCB/EERC-88/12Earthquake Engineering Research 
Center College of Engineering, University of California, Berkeley September 1988. 

[4] Hyo-Gyoung Kwak, Filip C. Filippou, "Finite Element Analysis of Reinforced Concrete 
Structures under Monotonic Loads”, Report No. UCB/SEMM-90/14, Structural Engineering, 
Mechanics and Materials Department of Civil Engineering University of California, Berkeley, 
Novemberl990. 

[5] Mander, J.B., Priestley, M.J.N., and Park, R. “Theoretical Stress-Strain Model for Confined 
Concrete", Journal of Structural Engineering”, American Society of Civil Engineers, Vol.114, 
No. 8, 1988, pp. 1804-1825. http://dx.doi.org/10.1061/(ASCE)0733-9445(1988)l 14:8(1804) 

[6] Sifat Sharmeen Muin, "A Parametric Study of RC Moment Resisting Frames at Joint Level 
by Investigating Moment -Curvature Relations”, International Journal of Civil and Structural 
Engineering, Vol. 2, No. 1, 2011. 

[7] Varghese P. C, “Advanced Reinforced Concrete Design ", PHI Learning Private limited. New 
Delhi, 2005. 


Moment Curvature 
Characteristics for 
Structural Elements 
of RC Building 


29 





